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The  Collinear  Antenna  Array. 
Theory  and  Measurements 
by 


Howard  W.  Andrews 
Ci  *ff.  ^ ^oratory,  Harvard  University 
Cambridge,  Mas sachusetia 


Abstract 


The  antenna  array  considered  herein  is  an  arrangement  of  a simple 
driven  dipole  and  two  parasitic  elements  with  axes  collinear  to  that  of  the 
driven  element.  A theoretical  solution  to  the  problem  is  obtained  by  con- 
sidering the  air  gap  between  the  elements  to  be  a simple  lumped  capaci- 
tancerof  a size  depending  upon  the  spacing  between  the  elements.  In  this 
way  the  array  is  considered  to  be  a re  actively -loaded  center-driven  dipole 
of  overall  length  equal  to  the  total  length  of  the  array.  The  reactances  are 
replaced  by  equivalent  generators  and  then,  through  superposition,  the  array 
is  separated  into  the  sum  of  two  dipoles,  one  driven  at  its  center  and  the 
other  by  two  off-center  generators. 

The  current  distribution  and  the  driving  point  impedance  are  computed 
for  a variety  of  length  as  a function  of  the  gap  spacing.  Measurements  are 
given  to  check  the  validity  of  the  theory  and  its  range  of  practical  application. 


I 

INTRODUCTION 

The  collinear  antenna  array  is  an  arrangements  of  driven  and  para- 
sitic elements  all  of  whose  axes  are  located  on  the  same  straight  line.  The 
parameters  needed  to  describe  the  array  are  themumber  an(l  spacing  of  the 
elements,  their  diameters  and  lengths,  and  the  nature  of  the  driving  voltages. 
The  coupling  between  the  elements  may  be  the  free- space  coupling,  a lumped 
reactance,  or  a two-wire  line.  It  is,  of  course,  necessary  to  limit  the  para- 
meters to  be  varied  since,  even  in  the  simplest  case  of  a three  - element  array, 
there  are  about  twelve  such  variables.  The  use  of  an  image  screen  reduces 
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the  variable  parameters  considerably  since  the  array  mu3t  then  be  symmetrical. 
The  configuration  actually  discussed  here  is  a three  element  array  with  a 
«lice  generator  in  the  center  element  at  the  plane  of  symmetry  of  the  array. 

The  elements  are  all  of  the  s-:.me  radius.  The  parameters  to  be  varied  are 
the  lengths  of  the  central  and  outer  elements  and  the  air  gap  sparing  between 
them.  The  quantities  measured  in  the  array  are  the  current  distribution  and 
the  driving-point  impedance  as  a function  of  the  three  variable  parameters. 

1-4 

This  array  has  been  previously  investigated  by  a number  of  people. 

Carter*  in  1932  included  it  in  a paper  considering  the  impedance  character- 
istics of  several  types  of  pairs  of  linear  radiators.  He  determined  the  self 
and  mutual  impedances  after  assuming  a sinusoidal  current  distribution  on 
radiators  of  length  equal  to  multiples  of  haif-w&veiengths.  The  expressions 
consisted  of  sums  of  sine  and  cosine  integrals  and  natural  logarithms.  The 
computed  results  were  fairly  good  since  he  had  restricted  himself  to  arrays 
in  which  the  even  current  distributions  were  resonant. 

Harrison^  in  1945  considered  an  array  of  two  identical  elements  driven 
by  identical  slice  generators  at  their  centers.  This  is  the  configuration  that 
results  from  a vertically  polarized  dipole  erected  over  a conducting  earth. 

His  procedure  was  entirely  different  from  that  of  Carter's  in  that  an  expres- 
sion is  derived  in  whit_h  the  current  is  the  unknown  quantity  satisfying  an  in- 
tegral equation,  thatis.the  current  distribution  is  defined  implicitly  by  the 
expression  and  must  be  obtained  from  it.  The  mathematical  technique  used 
is  that  of  Hallen^  for  solving  a similar  integral  equation  for  the  current  on  a 
simple  dipole.  The  integral  equation  is  obtained  by  first  writing  the  trigono- 
metric solution  of  the  simple  differential  equation  for  the  vector  potential  on 
the  antenna  surface  The  vector  potential  is  also  available  in  the  form  of  the 
Helmholtz  integi  al  of  the  current  distribution  on  the  antenna  Equating  this 
integral  to  the  solutions  of  the  differential  equation  results  in  an  integral 
equation  for  the  current  Two  of  these  are  obtained,  one  on  each  half  of  one 
of  the  identical  elements,  Hallen's  method  of  successive  approximations  is 
then  used  to  arrive  at  an  expression  for  the  current  distribution  for  both 
symmetrical  and  antisymmetrical  driving  generators  The  two  integral  equa- 
tions are  not  immediately  solvable.  The  boundary  conditions  are  introduced 
as  well  as  some  assumptions  ar.d  approximations  concerning  the  current. 
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Finally  explicit  expressions  for  the  currents  are  obtained,  in  the  form  of  a 
complicated  infinite  series.  The  computation  of  this  series  would  be  difficult 
and  Harrison  in  the  same  paper  presented  a more  approximate,  but  more 
tractable,  procedure  for  handling  the  integral  equation.  He  assumed  some 
additional  symmetries  in  the  current  distributions  that  led  to  computationally 
more  convenient  results. 

3 

R W.  P.  King  in  19  50  chose  to  consider  the  three  element  case  with 
only  the  center  element  driven  on  the  basis  that  this  array,  symmetrical 
about  the  generator,  is  the  only  practically  useful  arrangement.  His  procedure, 
was  similar  to  Hallen's  in  that  solutions  to  differential  equations  for  the  vector 
potential  were  equated  to  integrals  of  the  current.  A series  of  assumptions 
was  then  made  concerning  the  current  distribution  among  which  was  neglecting 
the  charging  current  at  the  ends  of  the  elements.  Considerable  use  was  made 
of  the  fact  that  the  vector  potential  at  a point  is  primarily  determined  by  the 
current  in  the  immediate  vicinity  of  that  point.  The  distribution  on  the  center 
element  was  obtained  by  driving  the  outer  elements  so  that  the  currents  at  the 
centers  of  all  elements  were  the  same;  this  was  done  for  symmetrical  and 
antisymmetrical  currents.  The  distribution  on  the  outer  elements  was  obtained 
by  considering  even  and  odd  distributions  as  well  as  symmetrical  and  anti- 
symmetrical  currents.  With  these  various  conditions,  a series  solution  to 
the  integral  equation  was  obtained  using  successive  approximate ons . From 
this  followed  expressions  fpr  the  self-'and  mutual  impedances  of  the  various 
elements . 

The  major  limitations  in  this  theory  are  two.  The  first  is  that  the  effect 
of  chargeable  end  surfaces  is  not  considered;  that  is,  the  model  best  repre- 
senting the  theory  is  one  that  has  ends  upon  which  no  charge  may  accumulate 
and  has,  as  chargeable  areas,  only  the  longitudinal  surfaces  of  the  conductor. 
The  second  limitation  is  that  no  large  odd  currents  should  exist  in  the  para- 
site, that  is  the  parasite  should  not  be  of  such  length  that  odd  currents  are 
resonant. 

At  the  beginning  of  the  present  study  a series  of  preliminary  measure- 
ments were  taken  to  check  the  impedances  computed  by  King.  The  driving 
point  impedance  of  several  combinations  of  length  was  measured  as  a function 
of  the  spacing  between  the  elements.  The  spacing  war  varied  between  zero-- 
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that  is,  actual  contact  between  the  elements  --  and  a large  fraction  of  a wave- 
length. For  this  case  of  contact  between  the  elements,  or  of  zero  gap  distance, 
the  array  degenerates  into  a simple  center  driven  dipole  of  length  equal  to  that 
of  the  array.  At  the  other  end  point  of  infinite  spacing  the  array  is  only  the 
driven  element  by  itself.  The  plot  of  the  driving -point  impedance  is  often  an 
arc  of  a circle  between  these  two  points.  It  may  also  be  a spiral. 

It  was  immediately  noticed  that  the  driving -point  impedance  changed 
very  rapidly  with  gap  spacing,  in  fact  the  complete  variation  in  impedance  often 
took  place  in  less  than  1/100  of  a wavelength  and  always  in  less  than  1 1 10  of 
a wavelength.  The  most  rapid,  and  also  the  greatest,  change  occurs  when  the 
gap  is  di  a high  current  point  for  the  dipole  that  results  when  the  elements  are 
in  actual  contact.  The  action  of  the  gap  is  to  reduce  very  rapidly  this  current 
at  the  gap  position  to  a comparatively  low  value.  This  large  change  with  gap 
size  of  the  magnitude  of  the  current  at  this  point  is  reflected  as  a correspond- 
ingly rapid  change  in  the  driving -point  impedance  as  well  as  a similarly  rapid 
change  in  the  magnitude  and  shape  of  the  current  distribution.  A comparatively 
slow  and  small  change  in  the  driving -point  impedance  occurs  when  the.  gap  is 
at  a low  current  point  of  the  dipole  resulting  from  direct  contact  of  the  elements. 
Then  the  current  is  already  comparatively  sma:l  and  reducing  it  to  zero  does 
not  have  a profound  effect  on  the  driving -point  impedance  or  on  the  current 
distribution.  For  this  case,  that  is  a minimum  in  the  current  distribution  at 
the  gap  point,  a spacing  of  as  much  as  1/10  of  a wavelength  is  necessary  foi. 
a complete  variation  in  the  driving -point  impedance. 

A comparison  of  these  experimental  results  with  the  zero -order  case 
computed  by  King  is  plotted  in  Fig.  18.  It  may  be  noted  that  the  shapes  of 
the  curves  are  essentially  similar  over  the  sections  for  the  larger  spacings. 

The  displacement  of  the  curves  is  no  doubt  due  to  the  fact  that  only  the  zerotVi- 
order  computation  has  been  carried  out.  The  agreement  and  trend  is  quite 
poor  for  the  small  (less  than  1/20X)  spacings.  This  disagreement  is  probably 
due  to  the  presence,  in  the  physically  existing  array,  of  end  surface  upon  which 
charge  may  accumulate  in  addition  to  the  chargeable  cylindrical  surfaces  of 
the  antenna  elements  The  existence  of  such  end  conditions  is  not  considered 
in  the  quasi- one -dimensional  theory  which  determines  the  surface  effects  using 
a line  current  distribution  at  the  center  of  the  antenna  cylinders  and  c Cju  s id  e r s 
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only  the  coupling  between  the  cylindrical  surfaces.  For  small  spacings  in 
the  actual  array,  the  chargeable  end  areas  are  sufficiently  close  that  there 
*s  an  appreciable  capacitative  coupling  between  them  as  well  as  between  the 
cylindrical  surfaces.  Hence  such  a configuration  requires  that  in  addition 
to  the  coupling  between  the  cylindrical  surfaces  adequately  treated  by  King, 
there  be  further  introduced  the  effects  due  to  the  end  coupling 

The  above  two  effects,  namely,  the  very  rapid  variation  of  driving- 
point  impedance  with  spacing  and  the  poor  agreement  of  the  King  theory  for 
small  spacing,  supposedly  due  to  end  capacitance,  leads  to  the  thought  tnar 
a theory  for  the  collinea.r  array  should  include  capacitative  effects  in  the 
region  of  the  gaps  as  an  essential  part  of  its  character.  For  close  spacings 
the  gap  is  actually  so  small  that  the  end  coupling  not  included  in  the  King 
theory  could  be  well  approximated  by  an  additional  near- zone  lumped  capaci- 
tance. This  would  then  permit  the  charging  current  of  the  lumped  capacitance 
to  be  treated  as  a displacement  current  across  the  gap.  It  would  attempt  to 
explain  the  poor  agreement  for  small  spacing  as  being  due  to  neglecting  the 
end  capacitances. 

With  the  premise  that  the  gap  is  a simple  lumped  capacitance,  the 
collinear  array  then  becomes  a simple  dipole  with  a capacitance  of  var- 
iable si*e  in  series  with  its  current  at  the  appropriate  point.  The  half 
length  of  the  dipole  is  equal  to  the  overall  length  of  the  array;  the  points 
of  insertion  of  the  lumped  capacitance  are  at  the  positions  of  the  gaps  in 
the  array.  This  configuration  will  be  assumed  to  represent  completely  the 
collinear  array.  See  Fig.  1 for  this  and  the  following  succession  of  events. 

The  use  of  a capacitance  in  the  array  introduces  the  possibility  of 
using  the  compensation  theorem  of  network  theory.  This  states  that  an 
impedance  in  which  a current  I is  flowing  may  be  replaced  by  a constant- 
voltage  generator  with  an  external  potential  difference  equal  to  -IZ  with- 
out changing  the  current  conditions  in  the  network.  Performing  such  a 
substitution  as  this  causes  the  array  to  become  a triply  driven  dipole.  The 
generator  at  its  plane  of  symmetry  still  exists  and  there  are,  in  addition, 
two  generators  to  replace  the  two  gap  capacitances.  The  array  is  otherwise 
a continuous  dipole  and  the  problem  could  be  solved  as  such.  The  generators 
would  be  introduced  as  the  energy  producing  boundary  conditions  on  the  scalar 
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potendal  in  an  approach  similar  to  Hallen's  iteration  technique  or  as  done  by 
Storer  using  a variational  principle.  The  linearity  of  Maxwell’s  equations, 
however,  allows  the  use  of  the  superposition  theorem  and  poses  the  possibility 
of  solving  the  problem  by  using  at  one  time  only  those  of  the  three  generators 
that  are  found  to  be  convenient.  The.  results  are  then  superimposed  to  give 
the  final  results. 

The  choice  of  solutions  found  to  be  most  convenient  is  to  break  the 

triply  driven  dipole  into  two  symmetrical  arrangements.  One  of  these 

is  the  doubly  driven  one  of  Fig.  2,  in  which  the  dipole  is  excited  by  two 

identical  generators  equidistant  from  the  center.  The  other  is  the  classic 

singly  driven  case  with  a slice  generator  at  the  dipole  certer  The  first  of 

6 7 

these  has  been  considered  by  Taylor  in  an  extension  of  S rer’s  variational 

ieennique  and  by  King^1  using  the  iteration  of  an  integral  equation.  The.  singly 

driven  case  is  a degenerate  form  oi  the  other  and  is  obtained  when  the  spacing 

between  the  two  generators  becomes  zero.  It  has  been  considered  by  Taylor^* 

7 

and  Storer  in  a variational  approach,  and  King  and  others  using  various 
techniques.  The  approach  used  in  this  thesis  will  be  that  of  the  variational 
principle  in  which  the  driving -point  impedance  is  an  extremized  function  of 
the  parameters  of  a trial  current  distribution.  Having  the  two  solutions,  they 
will  be  combined  in  such  phase  and  amplitude  that  the  doubly  driving  generators 
will  appear  to  be  a capacitative  reactance  of  the  proper  magnitude. 

In  conclusion,  the  purpose  of  this  research  is  to  develop  a theory  of 
the  coilinear  array  based  on  the  series  loading  of  a simple  dipole  by  lumped 
capacitances.  Computations  will  be  then  made  of  the  results  and  measurements 
taken  to  verify  them.  The  use  of  a two-wire  line  as  a coupling  between  driver 
and  parasite  will  be  investigated  experimentally.  This  coupling  is  not  simply 
represented  by  a lumped  reactance  and  its  corresponding  voltage  generator, 
but  rather,  because  of  its  ability  tc  support  an  unbalanced  mode,  a current 
generator  is  also  required. 
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THE  COLLNEAR  ARRAY  AS  THE  SUPERPOSITION  OF  A DOUBLY 
DRIVEN  DIPOLE  AND  A SINGLY  DRIVEN  DIPOLE 


I 


DOUBLY  DRIVEN  DIPOLE  SINGLY  DRIVEN  DIPOLE 


FIG.  2 MATHEMATICAL  MODELS  FOR  DIPOLES 
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Chapter  IJ 

A THEORY  OF  THE  ARRAY 

1.  The  Array  as  the  Superposition  of  Two  Separate  Dipoles 

The  closely  spaced  three -element  c.ollinear  array  will  now  be  investi- 
gated on  the  basis  that  it  is  well  approximated  by  a cylindrical  dipole  series 
loaded  by  a purely  reactive  impedance.  The  reactance  will  be  that  associated 
with  the  capacitive  coupling  between  the  ends  of  the  driven  and  parasitic  ele- 
ments. The  reasons  for  considering  this  approach  have  been  presented  in  the 
previous  chapter.  Only  the  center  element  will  be  driven  and  the  spacing 
between  the  center  element  and  the  two  parasites  will  be  the  same. 

It  will  first  be  necessary  to  discuss  in  more  detail  the  transition  from 
the  array  itself  to  its  representation  as  the  superposition  of  a doubly  and  a 
singly  driven  dipole.  The  series  of  steps  are  shown  graphically  in  Fig.  1. 

In  Fig  la  is  the  array  itself.  From  this  follows  (Fig.  lb)  the  u«e  of  a 
lumped  capacitance  as  the  complete  gap  effect,  then  the  substitution  of  an 
equivalent  generator  for  the  lumped  capacitance,  next  the  separation  (Fig.  Id 
and  Fig.  le)  into  the  two  dipoles,  and  finally  (Figs.  If  and  lg)  the  phasing  and 
amplitude  adjustment  of  the  generators  so  that  the  desired  equivalence  is 
obtained. 

The  voltage  generator  substituted  for  the  reactance  has  a potential 
equal  to  -IZ  where  I is  the  current  in  the  impedance  Z and  is  the  potential 
rise  of  the  generator  in  the  direction  of  the  current  flow  The  separation  cf 
the  triply  driven  dipole  into  one  that  is  singly  driven  and  another  that  is  doubJy 
driven  is  possible  because  of  the  linearity  of  the  applicable  field  equations 

Now  consider  the  phasing  and  amplitude  adjustment  of  the  doubly -driving 
generators  so  that  they  will,  when  superimposed  with  the  singly  driving  gen- 
erator, be  ninety  degrees  out  of  phase  with  the  total  current  through  them. 

The  symbolism  will  follow  that  used  in  Fig.  1.  Ic(z)  the  total  current  on 
the  colliuear  array;  its  value  at  the  gap  position  a = g is  (jh  g ) and,  as  is 

true  of  ail  currents  used  here,  I (-z)  = I (z).  The  current  distribution  on  the 

c c 

simple  dipole  of  height  h due  to  the  single  generator  V at  its  renter  (Fig.  le) 
is  I (z),  its  value  at  the  gap  is  I (+  g).  The  current  distribution  on  the  dipole 
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of  height  h due  to  the  doubly -driving  off-center  generators  (Fig.  Id)  at 

z = +p  is  I fz);  its  value  at  the  center  is  1,(0). 

— ” or  c. 

It  will  be  found  useful  later  to  use  current  distributions  that  have  been 

normalized  tc  a unity  driving  current  at  the  generator.  For  such  a unity 

cxirrent  at  the  center  of  the  singly  driven  dipole,  S(z)  is  the  resulting  current 

distribution.  Then  S(0)  = 1 and 


V 

(2-1) 

s 

where  Z is  the  driving -point  impedance  for  the  center  slice  generator  by 
s 

itself.  Similarly  for  the  off-center  generator  of  Fig.  Id  the  normalized 
current  distribution  is  denoted  by  D(z);  then  D(+g)  = 1 and 


Id(z)  = Zib(z> 

where  is  the  driving -point  impedance  presented  to  the  off-center  generators. 
The  desired  phase  and  amplitude  relation  at  z = _+g  when  the  two  are  super- 
imposed is 


- ic(±g)z 


V 


d 


or 


V±g) »- Vz 


(2-2) 


where  Z is  the  impedance  placed  in  series  with  the  current  on  the  loaded 
dipole  and,  in  this  problem,  ■will  be  the  apparent  reactance  of  the  gap.  At 
the  gap  position  (the  + and  - signs  will  be  omitted  from  this  point  on) 


y g)  --  yg)  + yg) 


is(g)  = Mg>  -yg) 

Substituting  (2-1)  and  (2-2)  above 


vs¥-  = - 


V,  v, 

d d 

~Z~  ' Z , 


= v - — — ) 

d’  Z Z,  ' 
d. 


1 
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Solving  for  Vj 


vd=va'¥"TTTl 


z z\ 


and  for  a unity  driving  voltage  V 


8 z +z 


Now  with  this  value  ol  the  total  current  lc(z)  at  any  point  on  the  array  is 


yz)  = is(z)  +id(z) 


= Ig( 0)S{z)  i Id(g)D(z) 


Vs  Vd 

■jr—  S(  z)  T ^ — D(a) 

s a 


/ 

= 2-S-S(z) 


for  V - 1/0° 
s — 


y*>  = 


V 

s 

2T 

B 


s<‘>  - 


Ic(-) 


D(z 

TTZ 


) 

1 


(2-3) 


The  driving-point  impedance  is  the  reciprocal  of  this  evaluated  at  the  driving 
point  of  z = 0,  thus 


Z 

c 


But  S(p)  = 1 so  that 


P(0)_S<g) 

1 + V * 


(2-4) 
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There  remains  now  the  evaluation  of  the  current  distribution  D{x) 
on  a dipole  of  height  b driven  by  identical  voltage  slice  generators  at  z = +g 
such  that  unity  current  flows  in  them.  The  arrangement  of  the  single  alice 
generator  at  the  center  will  be  obtained  by  setting  g equal  to  0 The  above 
algebra  will  be  used  to  superimpose  the  two  currents  in  order  to  obtain  the 
total  current  distribution  and  the  driving -point  impedance  of  the  array. 

2 . The  Stationary  Expression  for  the  Impedance 

The  technique  to  be  used  to  determine  the  driving '.point  impedance  and 

* 6 

the  current  distribution  will  employ  a variational  principle  for  the  impedance. 
From  this  will  be  obtained  the  current  distribution  by  means  of  the  Ritz  method 
applied  to  the  stationary  integral. 

In  general,  the  variational  method  depends  upon  obtaining  a functional 
expression  for  the  driving-point  impedance.  This  functional  expression  of 
a trial  current  distribution  is  so  formulated  that  its  value,  the  impedance,  is 
stationary  with  respect  to  small  deviations  of  the  trial  current  from  the  true 
current  distribution.  The  dependence  of  the  impedance  on  the  trial  current 
is,  in  fact,  only  of  second-order;  that  is,  errors  in  the  impedance  vary  as  the 
square  of  errors  in  the  assumed  current  distribution.  This  can  be  shown  by 
introducing  into  the  expression  for  the  impedance  a small  deviation  51  in  the 
current  distribution.  It  is  found  that  the  coefficient  of  the  resulting  error  6Z 
in  the  impedance  is  zero  bur.  that  the  coefficient  of  (5Z)^  is  not.  Alternatively, 
the  Euler -Lagrange  equations  for  the  expression  could  be  evaluated  and  it 
would  be  found  that  they  would  be  satisfied  by  known  expressions  containing  the 
Helmholtz  integral  for  the  vector  potential  and  the  definition  of  the  impedance. 

Once  a stationary  expression  for  the  impedance  has  been  found,  a 
trial  current  distribution  containing  a number  of  parameter  s io  ouuoi.ii.uted 
into  it.  Sinc“  the  expression  for  Z is  stationary,  the  Ritz  method  for  evaluating 
the  parameters  may  be  used  to  obtain  an  expression  for  the  current. 

Throughout  the  following  discussion  all  instantaneous  electromagnetic 

quantities  that  vary  in  time  are  understood  to  be  the  real  parts  of  complex 

quantities,  for  example.  H is  the  real  part  of  Ke^1.  The  time  dependence 
jwt 

e and  the  superscript  bar  denoting  a complex  quantity  are  omitted  throughout. 
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To  obtain  the  physically  existing  quantity  corresponding  to  any  equation  or 
expression,  multiply  through  by  e^Wt  and  take  the  real  part  of  the  result. 
Rationalized  MKS  units  will  be  used  throughout  this  discussion. 

The  mathematical  model  to  be  considered  is  that  drawn  in  Fig.  2a.  The 
antenna  is  assumed  to  possess  infinite  conductivity  and  to  have  only  currents 
in  the  z direction.  This  treatment  will  assume  no  currents  on  the  end  caps; 
this  condition  arises  from  the  approximate  one-dimensional  manner  in  which 
the  vector  potential  is  calculated  on  the  longitudinal  surfaces  of  the  antenna. 
Such  an  approach  is  justified  since  the  experimental  realization  is  an  antenna 
of  identical  axial  length  but  with  hemispherical  end  caps  of  the  same  total 
chargeable  area  as  the  corresponding  region  of  the  mathematical  model.  In 
the  experimental  model,  as  in  the  mathematical  one,  the  currents  at  z = + h 
are  zero.  If  the  desired  boundary  condition  is  approximated  experimentally 
in  this  way,  the  effective  height  can  differ  from  h by,  at  most,  a distance  of 
the  order  of  the  conductor  radius  and  probably  only  by  something  much  less 
than  this. 

The  boundary  conditions  resulting  from  the  application  of  Maxwell's 
equations  to  the  surface  between  a perfect  conductor  and  free  space  are 

n x H = -t  (2-5) 

n x E = 0 (2-6) 

The  first  equation  (25)  with  6 = -r  for  the  free  space  region  gives 

Hg  -tz(a,z) 

where  6 (a,z)  is  the  surface  current  density  in  the  z direction.  Note  that 
z 

a 2ir 

z t (r,0,z)rdrd9 
z 

= £ 2xa  6(a,z)  ( 2-7) 

X 

and  hence 

Id(z)  = 2?aHg(a,z)  -h^z^h 

The  second  boundary  condition  (2-6)  yields 
Ez(a,z)  = Ep(a,z)  = G 


(2-8) 
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Equation  (2-8)  does  not  hold  at  the  slice  generators  assumed,  to  drive  the 
antenna.  At  these  points  there  are  discontinuities  in  the  slope  of  the  scalar 
potentials.  These  slopes  will  be  assumed  to  give  a delta-function  discontinuity 
in  the  electric  field  such  that  in  the  limit  of  a discontinuity  of  infinitesimal  thick- 
ness 


A** 

/ Ez(a,z)dz 

J-g-t 


(2-9) 


Let  the  following  integral  definition  of  the  delta  function,  which  is  zero  every- 
where except  at  z = g and  there  takes  on  a value  such  that 


f(z)6(z-g)dz  = f(g)  , 


Tg- 


be  introduced.  Then 


'2-10) 


E (a,z)  = - V 6(z+g)  -hi£z<h  (2-11) 

Z & 

Note  that  (2-9)  has  E (a.z)  = 0 along  the  conductor  as  it  should  and  has 

z 

localized  the  generator  at  slices  at  a = + g. 


The  free-space  equation  for  the  electric  field  in  terms  of  the  scalar 
and  vector  potential  is 


E = - 


i>  . 


(2-12) 


Using  a periodic  time  dependence  and  considering  only  the  £- component  of 
E at  a distance  from  the  axis  equal  to  the  radius  a of  the  conductor  (2-12) 
become  s 

Ez(a,z)  =-jWAz(a.z)  - (2-13) 

The  vector  potential  A is  defined  by 

V x A = B 

the  Lorentz  condition 

V.  A = 0 , (2-14) 

c 


and  the  Helmholtz  integral 
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h 2tt 


A(a,9,z)  = £ 


f f 

^-h  ^6 


P (a.Q’.z')- adO’dz* , 

^ Z I — k.  -A.  i 


r -r 


(2-15) 


where 


|r-r'|  = ^(z-z'^+r^-fa*'  - 2ra  cos(9  - 9’)  , 


t (a.O1  ,z‘)  = z 6 (a, 6*  ,z') 

Z Z* 

For  the  case  considered  here,  where  the  antenna  is  cylindrically  symmetrical 

12-14 

and  of  sufficiently  small  radius  that  pa  « 1,  it  has  been  shown  that  a. 

good  approximatior\for  A,  even  at  the  radius  a and  z near  to  z1 , is  obtained 
by  assuming  all  of  the  current  to  be  located  at  the  axis  of  the  conductor.  Then, 
using  (2-7),  A at  the  radius  a is 


where 


Id{z')K(z.z')dz' 


K(  z,z! ) 


R = V ( z-z*  )^+a ^ 
= e / R . 


(2- 16) 


Note  what  has  taken  place  in  this  last  step.  In  (2-15)  the  vector  potential 
was  defined  exactly  as  the  integral  over  a tube  of  current  on  the  antenna  surface. 
In  (2-16)  this  was  replaced  by  an  approximation  using  a one -dimensional  current 
at  the  center  as  has  been  shown  reasonable  hy  King  and  Oseen.  As  a result, 
the  electric  field  (which  is  given  by  (2-lp)  and  (2-19)  below)  is  precisely  that 
at  a radius  a due  to  a current  at  r = 0. 

— A 

Since  A has  only  a z component  (as  seen  from  the  integral)  then 


^ 9A  (a,z) 

V • A(a,z)  = 

'd  z 

Substituting  this  in  the  Lorentz  condition  (2-14)  gives 
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3A  (a,z) 

+ ^(a,z)  = 0 

a z c 


and  hence. 


,,  v • 9 A (a,z) 

8p(a,z)  _ joj  zv 

3 z c 8 z2 


Then  (2-12)  tor  the  electric  field  at  the  radius  a becomes 

a2A  (a,z) 

Ez(a,z)  = - juA^a.z)  - — z 


a z 


(2-17) 


where  |3  = w/ c . Rearranging  (2-17)  leads  to 


Ez(a,z)  = + pZ)A_(a,z) 

£ 8z 


(2-18) 


= L(z)Az(a,z)  (2-19) 

By  equating  (2-19)  for  the  electric  field  at  a distance  a from  the  one 
dimensional  current  Id(z)  and  (2-11)  for  the  field  on  the  surface  of  a perfect 
conductor,  there  is  obtained 


- Vd6(z  + g)  = L(z)Az(a,z)  (2  20) 

Multiplying  (2-20)  through  by  l^(z)  and  integrating  from  -h  to  h gives  an 
equality  relating  the  complex  power  supplied  by  the  generator  to  the  complex 
power  radiated  by  the  current  z)  as  evaluated  at  r = a. 


h 

Id(z)vd6'z  i g)dz  = f 

^h 


Id(z)L,(z)Az(a,z)dz 


Performing  the  integration  on  the  left  using  (2-10)  and  substituting  (2-16) 
on  the  right  give  8 


h 

r 


-y-^v^d^  / dzyz>Mz) 


I^z'JKtz.z^dz' 


(2  21) 


Substituting  the  definition  of  the  impedance. 
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Vd_Id'ig)Zd  * 

into  (2-21)  and  making  vise  of  the  symmetry  of  the  current,  I (-g)  = I (g>. 

z * z 

leads  to 


h h 

- J dzId<">L(Z>  f 

-h  -h 


I^(z')K(z.z')dz' 


(2-22) 


The  presence  of  a minus  sign  on  the  left  locates  the  power  source  in  the 
generator;  the  right  side  is  the  power  radiated.  The  solution  of  (2-22)  for 


Z,  is. 
d 


d Siring) 


h h 

j£ P dz  Id(z)L,(z)  > Id(z')K{z,z')dz' 

4 4 


< 2-23) 


This  expression  for  Zd  is  stationary;  that  is,  small  deviations  of 
the  values  of  Ij(2)  from  the  true  value  will  result  in  only  second-order 
variations  in  Z^.  This  can  be  shown  by  obtaining  the  first  variation  6Zd 
of  Zd  due  to  a small  variation  61  in  Id(z).  The  volue  of  SZd  is  found  to  be 
zero. 

The  expression  (2-23)  can  be  simplified  using  the  symmetry  of  Id(z) 
and  L(z)A  (a,z).  Using  these  properties  (2-23)  becomes 


Zd=- 


J± 


b 

J dzId(z)L(z)  JT  Id(z')K(z,z')dz'  (2-24) 


4 


4,rId(g)  J0 

In  (2-24)  L(z)  operates  on  the  kernel  of  tlie  integral  Since 

2 r.2 

d K{z,zl)  = -^-j-Ktz.z')  , 

8z  8z' L 


Id(h)  - 0 = Id(-h)  , 

and  ^d(z)  and  L.(z)Id(z)  are  even  functions  about  z - 0,  it  is  possible  to  simplify 
the  L(z)A(z)  part  of  (2-24).  The  result  is 
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f-  h h 

f dzl,(z) 

1 d 


h 

i [K(z,z>] 

'•'O 


+ K(z,- z’l]  L.fz’H^z’Mz' 


dl  <z) 

(3^  dz 


) / Id5a)[K(a,h)  + K(z,~b.];]dz 

Z = h v/, 


(2  25) 


3 The  Trial  Current  Distribution 

The  selection  of  a trial  function  to  be  used  in  (Z- 25)  is  governed  as  much 
by  the  necessity  of  obtaining  results  in  terms  of  tabulated  functions  as  by 
employing  an  excellent  approximation  to  the  c urrent  distribution.  With  a 
kernel  of  the  type  occurring  in  (2-25),  trial  functions  including  sinpz,  cospz, 
z sinpz,  zcospz  and  constants  yield  the  tabulated  generalised  sine  and  cosine 
integrals.  The  distribution  is  known  to  be  continuous,  even  about  z - 0 zero 

at  r h:  and  to  have  a discontinuous  slope  at  z =+  g.  All  of  these  conditions 

““  ' 7 

should  be  approximated  as  closely  as  possible.  Storer  found  that  a combination 

of  3inpz,  cospz,  and  a constant  quite  accurately  represented  the.  simple  dipole 

current.  Such  functions  should  also  be  suitable  for  section  of  the  antenna 

defined  by  |zj3>g.  A constant  and  cospz  should  be  satisfactory  for  !z|<.g. 

With  these  conditions  ir*mind  suitable  trial  currents  are. 

I (z)  = C,  + C,  cosBz  |z|<g 

Q.  1 L.  

y*)  - k>1  sin  p(h-  z)  + D2[  1 - cosp(h-z)]  \z\ > g 

where  the  C’s  and  D!s  are  complex,  coefficients.  Note  that  Z in  (2-23)  is 
independent  of.  the  absolute  value  of  the  current  levels  Hence,  letl^l+g)  = 1. 

The  trial  currents  then  become 

ld(z)  = 1 + C [a  + cospz]  |z|  <g  (2-26) 

I^Cz)  = 6[  1 - cosp(h-z)]  + D^sinp(h-z)  + e [1  - cosp(h-z)]j  |z|>g 

(2-27) 


where 
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Q = - COS  Pg 


1 - cos  p(h-g) 

£ = ..  sin3(h  a) 

1 - cos  p(h-g) 

The  lengths  for  which  these  may  be  considered  reasonable  are  pg  <Cir 
for  (2-26)  and  p(h-g)  < 2ir  for  (2-27),  Equation  (2-26)  with  a constant  term 
would  be  expected  to  fail  for  g near  a half- wavelength  Equation  (2-2?j 
actually  has  a singularity  at  (h-g)  equal  to  a wavelength  and  probably  is  not 
too  accurate  for  (h-g)  larger  than  about  thr ee -quarters  of  a wavelength. 


4 . Evaluation  of  the  Impedance  and  the  Current  Distribution 
Substituting  (2-26)  and  (2-27)  into  (2-25)  yields^ 

zd=4T^o  + Ycc  +ydd  + Yccc2  + Y DDD^  +yCDCD^ 

(2-28) 


where 


[ l-cos  p(h-g)] 


[ 1 -f  cos2p(h-g)j  sinpa-cos2p(h-g)sinpy(2g)Z+a' 


- 2 cos  p(h-g)[  sin  $\f  (h-g)2+a2  - sinp^h+g)  2+a2]  - sin  py(2h)2+a2 
+ f 2Pg  c.os2p(h-g)  + sin  2p(h-g)j  ’C(pa,2pg) 


I [ 2p(h-g)cosp(h-g)  - 2 sinp(h-g)]  C(pa,p(h-g))  + 2PhC(pa,2ph) 

- [ ZpthTgicos  p(h-g)  x 2 sinp(h-g)]  C(pa,p(h+g))  + 2Cs(pa;p(h-g)) 
+ 9in2ph[2Cc(pa,p(h+g))  - Cc(pa,2ph)  - Cc(pa,2pg)] 
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- cos  2ph[  2Cg(patp(h+g))  - cypa.zph)  - Ca(pa.2pg)] 

- j(coe2p(h-g)  cob  p^/  (2g/ba  - [ 1 + cos“p(h-g)]  cospa 

+ 2 coc  p(h-g)[  cos  $\J  (h-g)2+a2  - cos  $\J  (h+g)2+a2]  + cos  $\J{2h)^  +2.^ 
+ [ 2p  cos2  p(h-g)  + sin  2{3(h-g)’]  S(pa,2pg) 

+ [ 2p(h-g)cos  B(h-g)  - 2 sinp(h-g)  S(j5a,p(h-g))  i 2phS(pa,2ph) 

- [ 2p(h+g)cos  p(h-g)  + 2 sin  p{h-g)]  S(pa,p(h+g))  + 2Sg(pa,p(h-g)) 

+ sin  2Sji[  2S  (pa,p(h+g)>)-  S (pa,2ph)  - S (pa,2pg)] 

C ^ 

*1 

- cos  2ph[  2Sg(pa,p(b+g))  - Sg(pa,2ph)  - Sg(pa,2pg)]  • , 

= - , coapgcos  p(h-g)[  sin  pa  - G:npV(?g2)  + a2  ] 

1 -cos  p(h-g) 

- cos  pg[  sin. p^/  (h-g)2+a2  - ainp\J  (h+g)2+a2j 

+ [ 2p  g cos  pg  cos  P(n-g)  + sin p (h-2g)]  C(pa,2pg) 

+ [sinpg  + p(h-g)cospg]  C(pa,p(h-g)) 

+ [sinpg  - p(h-fg)  cos  pg]  C(pa;ft(h+g)) 

+ sinpb["Cc  (pa,p(h+g)>  - ’Cc(pa,p{h-g))  - Cc(pa,2pg)] 

- cos  ph[  Cg(pa,p(h-fg))  - Cs(pa,p(h-g))  - Cs(Pa,26g)] 

- j{cos  Pg  cos  p(h-g)[  cccp^  (2g)  -t-a  - cospa] 

- coS?g[cosp\^V  - cos  p \j  (b-g)2+a2] 
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+ [ 2pg  cos  pg  cos  p(h-g)  + sin  p(h  - 2g)]  S(pa,2ps) 

+ [sinpg  + p(h-g)cospg]  S(pa.p(h-g)) 

+ [ sin  Pg  - p(h+g)  cos  pg]  S(PaTp{h+g>) 

+ sinph[Sc(pa,p(h+g))  - Sc(pa,p(h-g))  - Sc(pa,2pg)l 

- cos  ph[  Sg(pa,p(h+g))  - Sg(pa.p(h-g))  - Sg(pa,2pg)j)  j , 

J 

» 

V = y 4 sin  2p(h-g) sin  p\/ (2g)2+a2  + sin p(h-g)sin py(2h)2+a2 

D [ 1-cos  R(h-g)]  2 [2 

+ [ sinp(h-g)  + ~ in  2p{h-g)]  [ sinpV^h-g)2+a2  - sin p^(h+g)2+a2  - sinpa] 

- [ 2ph  sinp(h-g)  + cos  p(h-g)  - l]  Qpa,2ph)  - 2 sin p(h-g)Cg(pa,p(h-g)) 

- [ Pg  sin  2p(h-g)  - cos  2p(h-g)  + cos  p(h-g)]  C(pa,2pg) 

+ [ P(h+g)(sinp(h-g)  + -^sin  2p(h-g)H2  * 2cosp(h-g)]  C(pa,p(h+g)) 

- [ P(h*g)(sinp(h-g)  + j sin2p(h-g))  - 2sin2p{h-g)]  CCpa.pOi-g)) 

[cos  p(h+g)  - cos  2ph]  [ 2Cc(pa,p(h+g))  - Cc(pa,2ph)  - TTjpa.Zpg)] 

- [ sin  p(h-fg)  - sin2ph][2Cg(pa,p(h+g))  - Cg{pa,2ph)  - Cg(Pa.2pg)] 

- j(-  sin  2p(h-g)cos  p\/V2g)  2+a2  - sin  p(h-g)  cos  P V(2h)2  +a2 

- [ sinpTh-g)  + -^  sin  2p(h-g)]  [ cos  p^(h-g)2+a2  - cos  p \J  (h+g)2+a2  - cospa] 

- [ 2ph  sin  p(h-g)  + cosp(h-g)  - l]S(pa,2ph)  - 2 sin  p(h-g)Sg(pa,p(h- gl) 

- [pg  sin  2p(h-g)  - cos  2p(h-g)  + cos  p(h-g)]Sipa,2pg) 
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t f p(h-f-g)(sinp(h-g)  + sin  2p(h-g))  + 2 - 2cosp(h-g)]  S(pa,p(h+g)) 

- [ P(h-g)(sinp(h-g)  + -|  sin  2p(h-g))  - 2 sin**  p(h-g)]  S(pa,p(h-g)) 

- [ cos  (h+g)  - cos  2£h]  [ 2Sc(pa,p(h+g))  - Sc(pa,2ph)  - S^pa^g)] 

- [ sin  p(h+g)  - sin2ph]  [ 2S  (pa,p(h+g))  - S (pa,2ph)  - S (pa,2pg)])  • 

a o o 

Y CC  ~ cos2Ps[  sin  pa  - sinpV(2g)Z  + a2]  + Cs(pa.2pg) 

+ 2 cos pg[  pg cos  pg  - sinpg]  C(pa,2ps) 

- j(  cos2  pg[  cos  py (2g)Z+aZ  - cos  Pal  + S^pa.Zpg) 

+ 2cos  pg[  Pg  cos  pg  - sin Pg]  S(p'a,2pg))  , 


DD 


1 - cos  p(h-g) 


(1  + cosp(h-g))[  sin  pa  - sin  p,/(h-g)Z+a?' 


+ sin p/|/(h+g)Z+aZ  - ^ sinp^/(2h)ZraZ 


- j sinpy (2g)Z+aZ]  + 2Cs(pa,p(h-g)) 


+ [ 6g(  1 + cos  p(h-g))  + sinp(h-g)]  C(pa,2pg) 


+ [ph{l  + c.osp(h-g))  - sinp(h-g)]  C(pa,Zph) 


+ [P(h-g)(l  + cos  p(h-g))  - 2 sin  p(h-g)]  C(pa,p(h-g) 


- P(b+g)(i  + cos  p(h-g))  C(pa,p(n+g)) 


+ sinp(h+g)[2Cc(pa,p(h+g))  - ^(Ba.Zph)  - ^c(pa,2pg)] 
- cosp(h+g)[  2C  (pa.B(h-t-g))  - C (pa,2ph)  - C (pa,2pg}] 

O 5 b 
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- j((l  T cos  P(h-g))[  COB  py  (h-g)2*a2  - cospa  - cos  py (h+g)2+a2 
+ icospy(2h)2+a"2+  ^cospy(2g)2+a2]  + 2Ss<pa>p(h-g» 

t [pgtl+cos  p«h-g)H  sinp(b-g)]  S(pa,2Sg) 

+ [ «h{  1 + cos  p(h-g))  - sinp(h-g)]  S{pa,2ph) 

+ [P,h-g)U  f cos  p(h-g))  - 2 sinp(h-g)]  S(pa,p(h-g)) 

- P(h+g)(l  + cos  p(h-g))  S(pa.p(h+g)) 

+ sinp(h+g)[2Sc(pa,p(h+g))  - SJpa,2phl  - SJpa^pg)] 

' 

- cos  p(h+g)[  2Sg(pa,p(h+g))  - Sg(pa,2pn)  - Sg(pa,2pg>]  ) 

and 

Vnr  1 cos  Pg  sin  p(h-g)[  sin  pi / (2g)2+a2  - sinpa 

'-'u  1 - cos  p(h-g)  ^ 

+ sin  p|/<,h-gl2+at'  - sin  py^h+g)  2+a2] 

- [ 2Pg  cos  fig  sinp(h-g)  + cos  pg  - cosp(h-2g)]  C(pa,2pg)l 
+ [ P(H+g)  cos  pg  sin  p(h-g)  + cos  ph  - cos  Pg]  <5(pa,p(h+g)» 

- [pfh-glcospg  sin p(h-g)  - cos  Pg  + cosp{h-2g)]  C(pa,p{h-g)) 


+ [ cos  ph  - cos  Pg]  [ C (pa,p(h+g))  - C (pa.p(h-g))  - C fpa,2pg)] 


[ sin  ph  - sinpg][Cg{pa,p(hTg)>  - Cs{pa,p{h-g)}  - C^pa^pg)] 


- i { cos  pg  sinp(h-g)f  cos  pa  - cos  p f I 2g)2+a2  + cos  p \j (h+gl2+a 


os 


pl/(h-g)2  r a2] 
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-|  2£gcospg  sin p(h-g)  + cospg  - cosp(h-2g}]  S(Ba,2pg) 

+ [ p(h+g)  cos  Pg  sin  P(h-g)  + cos  Bh  - cos  Pg]  S(Ba,p{h+g)) 

- [ p(h-g)  cos  pg  sinp(h-g)  - cos  pg  + cosp{h-2g}]  S(pa,p(h-g)) 

+ [ cos  ph  - COS  pg][  Sc(pa,p{h+g))  - Sc(pa,p(h-g))  - Sc(pa,2pg)] 

+ [sin  ph  - sinpg][  Sg{pa,p{h+g))  - 2s(pa,p(h-g))  - Sg(pa,2pg)]  )|  . 

See  Appendix  A for  the  evaluation  of  the  various 

This  expression  (2-28)  is  known  to  be  stationary;  that  is,  it  is  an  extre- 
mum of  the  function  I^(z).  More  exactly  8Z/8I  = 0 or 

az  az  _ 
ac  " 8D 

After  carrying  out  these  operations  and  solving  the  resulting  simultaneous 
equations  for  C and  D the  following  results  are  obtained: 

c - yd7cd  "^cYdd 
4yCCVDD  ~ 1 CD 

D = YCVCP 

4vCC VDD  ' YCD 


5.  The  Single-Driven.D  pole 

The  singly-driven  dipole  is  a degenerate  case  of  the  doubly-driven 

one;  see  Fig.  2’-2b.  Allowing  g to  equal  zero  causes  the  two  slice  gener- 

* 6 

atom  to  become  hut  one  at  z=  0.  This  causes  v _ = v = v = 0 and 

' i * 'i.i.  • : n 


___1 

(1-  COSph)‘ 


einpa  - sin  p]/(2h)  ^+a^  - 4sinph  C(pa,ph) 


+ 2ph  C(pa,2ph)  + sin  Z0h[  2C(,(pa,ph)  1 2C„(pa, 2ph)] 


o 
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- cos  2ph[  2C  (pa.ph)  - C_(pa,2ph)]  + 2G  (pa,2ph) 

S o'  S 

- j(cos  py  ( 2h)^  + - cos  ^a  - 4 sin  ph  S<pa,0h) 

+ 2(3 h S(pa,2(3h)  + sin  2ph[  2Sc(pa,ph)  - Sc(pa,2ph)] 

- cos  2ph[  2Ss(pa,ph)  - Sg(pa,2ph)]  + 2Sg(pa,2ph))|  , 


^ ( 1 -cos  ph)^ 


s^n  ph[  sin  p y(2h)t'+a^  - sinpa]  - 2sinph  C (pa.ph) 


- 2, 


- [2phsinph  + cosph  - l]  C(pa,2ph)  + 2[  1 + sin  ph  - cosph]  C(0a,ph) 


[cosph  - cos  2ph][  2Cc(pa,ph)  - C^(pa,2ph)]  - [sinph  - sin  2ph] 


[ 2Cs{pa,ph)  - Cg(pa,2ph)] 


j(sin  ph[  cos  pa  - cos  p|/(2h)2  + a2]-  2 sinph  S (pa,ph) 


• 2, 


- [ 2ph  sinph  + cosph  - 1]  S(pa,2ph)  + 2[  1 + sin  ph  - cosph]  S(pa,ph) 


- [cosph  - cos  2ph]  [ 2$  (pa.ph)  - Sc(pa,2ph)]  - [sinph  - sin  2ph] 


[ 2S  (pa.ph)  - S(pa,2ph)]} 


and 


DD 


4(1  + cos  ph)[  sinpa  - sin pl/( 2h)^+a - 2 sinph  C(pa,ph) 

l-cosph 


+ 2Cg(pa,ph). 

+ [ ph(  I + cosph)  - sinph]  C(pa,2ph)  + sin  ph[  2Cc(pa,ph)  - CT(pa,2ph}] 
- cos  ph[  ZC^Cpa.ph)  - C (pa, 2ph)] 

■j  3 
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Chapter  III 

METHOD  OF  IMPEDANCE  MEASUREMENTS 
AND  AU  X3  L LI  ARY  MEASUREMENTS 


1 . Line  Theory 

The  impedances  presented  in  this  paper  were  measured  on  the  coaxial 
transmission  line  with  a characteristic  impedance  equal  to  125  6 ohms, 
a phase  constant  p of  12.  775  radians  per  meter,  and  a theoretical  attenua- 
tion constant  a of  0.003  nepers /mete r. 


The  differential  equations  describing  the  current  and  voltage  at  points 
remote  from  the  ends  of  the  line  are 


ItJ—L  = (r  + jw&)I(z) 


(3-1) 


-ai(z) 

3z 


fg  + jux)V(z) 


(3-2) 


where  r,6,g  and  c.  are  the  series  resistance  and  inductance,  shunt  cor.  due  - 
tance  and  capacitance  per  meter  cf  line.  The  lime  dependence  and  tie 
complex  nature  of  the  current  and  voltage  have  been  suppressed.  A soiu 
tion  for  the  current  containing  the  boundary  conditions  indicated  in  the 
diagram  above  is,^ 

r 

1(1)  = V !- 


Z_  sinhy  w + Z^  coshy  w 


T~~ 


^ ( Z(  + Z^Z^sinhys  + Z{  ( Z^  + Zg)cosh  y : 
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wher  e 

« y Ti 

■y  = a + jp  = [ (r+jujtHg+jwc.)]  A/ 41 

Using  the  form 

Z = Z^coth  9 = Z^coth(p+j$)  , (3- 3) 

an  expression  containing  only  hyperbolic  functions  is  obtained. 

.r  ["  sinhO  sinhlyw+O  ) 

'Z'  “ Zc  j "smh(ys  + 0o+0s) 

7.'his  formula  describes  the  current  completely,  but  is  in  complex 
form.  A detector,  sensitive  only  to  the  magnitude  of  this  current,  will 
measure  the  quantity. 


«■>!  * £ 


r 2 2 2 2 

j (sinh  pQ4-sin  ^Q){sinh  (a,u>+p  g)  + sm  'Pw+6^)} 


2 ” 2 

sinh  (us+p^fp^l  + sin 


1/2 


The  position  of  the  minimum  of  this  current  distribution  may  be  conveniently 
located,  and  for  a lossless  line  this  position  is  solely  a function  of  the  phase 
4>  of  the  load.  Taking  the  derivative  of  the  square  of  the  magnitude  and 
equating  to  zero  results  in  the  equation 


a sinh(aw+p  )cosh(aw+p  ) = - p sin(pw  + J>  )cc3(pw+^  ) 
s s s s 


Squaring  both  sides  and  using  the  double  angle  trigonometric  identities 
gives. 


and  therefore, 


As  suming 


a2  2 2 

sinh"2law+p  ) = 1 - cos  2^pw+i6 

S £ 


A / ^ 

C 08  2(pw+^  ) = +1  / 1 - sinh^2(iaw+p  ) 

8 “!/  S 


“ -y 

— ^ sinh  2(  uw+p  ) 1 

P 


- 


<3-4> 
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then 


cos  2((Jw  + ^g)  = + 1 


and  hence 

pw  + 4>&  - jh  ~ , n = 0,1, 2, 3,.  . . . 

determines  the  positions  of  the  maxima  and  miaima.  The  minima  occur 
for  n even.  For  the  first  minimum  from  the  end  (w  = 0),  n = 2,  w = wmin 
of  the  current  distribution  and 


i 


s 


TT  - pW 


min 


(3-5) 


The  other  parameter  characterizing  the  load  impedance  is  pg.  It 

may  be  determined  from  the  current  distribution  in  the  vinicity  of  the 

minimum  normalized  to  the  value  of  this  minimum  current.  There  are  two 

points,  one  either  side  of  the  minimum,  having  the  same  value  of  current. 

2 1 

For  a lossless  line  ((a6w)  <V— ) these  points  are  located  svmmetrically 

O 

at  a distance  6w  either  side  of  the  minimum.  In  the  vicinity  of  the  minimum 


u&y 


min' 


= P(w) 


sinh{yw  + p q6w) 

sinh{yw  + 9 ) 

' T m s' 


The  squared  magnitude  of  this  quantity  is. 


2 2 

, sinh  (aw  + p + a&w)  + sin  (fjw  + 4>  + (35w) 

p ( w)  


m r s 

~~T2 


inh  (awm  + p a)  + sin  (pwm  + 6J 


m 


Substituting  in  the  value  pw  + ^ = it  at  the  first  minimum, 

8 

p*”(w)  sinh  (aw  + p ) = einh  (aw  + p + a.5  ) + ?in  p5v.r 

Hi  6 m s — ni 


2 2 2 

Assuming  that  a6w«|awm  + pgI  then  >p“(w)  - 1)  sinh  (aw^  + og)  = sin  p6w. 

The  point  of  maximum  slope  of  p2(w)  occurs  at  p2  = 2,  and  it  is  the  point 
where  the  current  squared  on  the  line  is  twice  that  at  the  minimum.  For 
p2(w)  = 2 

p = sinh  *(sin6&w)  -aw  (3-b) 

s m 

Thus  the  two  quantities  needed  in  (3-3)  to  describe  an  arbitrary  load  impedance 
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4.  xv  x t O 

ran  be  determined  from  the  position  of  the  minimum  and  the  width 
of  the  distribution  curve  at  the  double -power  points. 

Three  approximations  have  been  made  here,  other  than  those  in 

(3-1)  and  (3-2),  in  locating  the  minimum;  (3-5)  and  \y  = )/P 

are  true  to  the  extent  that  (3-4)  is  satisfied.  For  the  experimental 

line  used,  and  mr  example  with  an  unusually  lew  standing  - wave  ratio 

of  1,1,  then  {—^)  sinh^2(aw  + jj  ) = 5.  54  x 10  ^ and  (1*  4)  is  certainly 

well  satisfied.  Using  the  formula  Q = p/2a  = 4580,  the  condition 

«o6w)^<^  i becomes  Q^^>^>1,  and  a5w  <«  aw  + p requires  Q>'j>l.  All 
o m s 

of.  the  conditions  on  the  approximations  are  easily  met  ?,nd  the  line 
losses  do  not  contribute  to  the  error  in  (3-  5)  or  (3-6). 

The  limiting  accuracy  thus  rests  on  the  accuracy  with  which  the 

2 

position  of  the  double  power  points  may  be  located.  The  choice  of  p =2 
vs  optimum  for  a square  law  detector  in  that  the  points  for  w^  6w  fall 
on  the  steepest  part  of  the  measured  distribution.  Attempts  to  average 
a few  points  about  the  minimum  as  compared  to  averaging  the  position 
of  the  two  double  power  points  resulted  in  no  improvement  in  the  accuracy 
'with  which  the  minimum  was  located.  The  error  in  measuring  the  double 
power  width  increases  at  lower  standing-wave  ratios,  but  will  be  less 
important  since  the  impedance  corresponding  to  a small  standing -wave 
ratio  is  then  less  sensitive  to  errors  in  the  determination,  of  p . Sinr.r 
only  a 6 db  range  of  the  detector  calibration  curve  is  used,  it  is  possible 
to  restrict  the  measurement  to  that  portion  found  to  have  a constant  slope 
and  be  nearly  smia'-e  law. 

2 Computation  of  Impedances 

A desk  calculator  used  in  conjunction  with  a large  scale  Smith 
chart  provides  a convenient  technique  for  processing  the  data  from  the 
slotted  line  measurements  T.  sum  a,  u average  of  the  positions  of  the 
two  double-power  points  are  quickly  computed  to  give  the  positions  of.  the 
minimum,  and  thus  6 . Their  difference  gives  Sw  and  thus  p . Tables 
of  svnh  n,sinp6w)  for  0 <C  p6w<^  1.  500  in  steps  of  0.000  5 radians  were 
prepared  to  facilitate  the  computation  of  pg.  The  correction  due  to 
c.w  in  (3-6)  was  of  the  order  of  0.001  radians  and  usually  negligible. 
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A 15  inch  diameter  Smith  chart  was  covered  with  a rotatable  transparent 

d sk  c arrying  a radial  p scale.  Knowing  p and  i>  the  normalized 

8 8 8 

impedance  on  the  line  was  calculated  with  about  1 percent  Accuracy . 

17 

Some  of  the  data  were  corrected  for  the  end  effect  on  a coaxial 
line.  The  admittance  corresponding  to  the  measured  impedance  was 
obtained  graphically  from  the  Smith  chart  by  rotating  the  measured 
p <6g  point  by  180  degrees  about  the  center  of  the  cha.-rt.  The  end - 
correction,  a negative  capacitive  susceptance,  of  0.  144  mhos,  normalized 
at  600  MC  on  th' s line  is  added  and  the  new  p and  ^ point  rotated  another 
180  degrees  to  give  the  corrected  measured  load  impedance.  See  Fig.  1 
and  Figs.  15  to  18  for  examples  of  the  importance  of  this  correction. 

• ^ g 

A more  complicated  correction  has  been  considered  by  Whinnery  and 
19 

also  by  Zeoli,  but  it  was  decided  that  the  general  accuracy  of  the  problem 
is  not  great  enough  to  make  such  correction  necessary.  The  total  capacity 
of  the’r  tt  network  is  nearly  equal  to  that  of  the  Hartig  correction  if  the. 
inductance  is  neglected. 

3 ..  General  Check  on  System  using  Half  Dipole  Impedance  Measurements . 

The  general  accuracy  of  the  measurements  was  checked  by  measuring 

the  impedance  of  the  simple  half  dipole  of  constant  radius  as  a function 

of  length,  correcting  it  for  end  effects  and  then  c omparing  the  results 

with  the  King  Middleton  second-order  impedance  for  a constant  radius 

20 

dipole,  obtained  recently  by  cross  plotting  from  constant.!!  data.  See 
Fig.  3 for  a comparison  of  these  results. 


The  accuracy  of  the  auxiliary  probe  as  used  on  the  polyfoam 
supporting  colum  was  investigated  by  measuring  the  current  distribu- 
tion on  a.  half  wavelength  dipole  with  both  probes.  The  resulting  curves 
were  plotted  and  normalized  to  have  the  same  max’ mum.  (See  F’g.  4 
Also  plotted  m this  diagram  are  the  King  Middleton^*  first-order  dis- 
tribution and  th*1  current  obtained  using  the  equation  ^or  I (z.)  on  page  2 24. 

s 
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5.  Measurement  of  Gap  Capacitance. 

Experimentally,  the  gap  between  the  end  of  the  driven  antenna 
and  its  parasite  is  the  parameter  -varied  in  any  particular  array.  A 
varying  capacitive  reuctcuiue  is  assumed  to  be  the  equivalent  of  this  in 
the  theory  of  the  array.  Hence  the  magnitude  of  the  lumped  capacitance 
between  the  hemispherical  end  caps  is  needed  in  order  to  compare  the 
theory  with  data  using  the  gap  reactance  as  a parameter. 

Since  the  theory  applies  precisely  only  for  small  (0  <1  gap  X / 1 0) 

gaps,  it  seems  reasonable  that  the  capacitance  of  such  a gap  could  be 

measured  by  constructing  it  in  the  center  conductor  of  a shorted  coaxial 

line.  Such  a measurement  was  set  up  with  a variable  gap  placed  a half 

wavelength  from  the  shorted  end  of  the  line.  The  apparent  series 

reactance  and  capacitance  of  thegap  was  so  determined  and  plotted  in 

Fig.  5 . For  spacings  less  than  0.  01  cm  the  curve  was  extrapolated 

linearly  from  its  measured  values  at  larger  spacings.  The  results 

22 

agree  qualitatively  with  those  found  by  Jeans  for  the  capacitance 
between  two  identical  spheres.  For  very  small  gaps  in  the  actual  array, 
a small  piece  of  dielectric  tape  was  used  to  separate  the  ends  of  the 
antenna  elements.  The  tape  was  0.006  cm  thick,  and,  assuming  a dielectric 
constant  of  2 or  3,  the  equivalent  air  gap  is  about  0.003  cm  with  a 
corresponding  gap  reactance  of  approximately  250  ohms. 


n 


A A MEASURED  IMPEDANCE 

• — • MEASURED  IMPEDANCE 
CORRECTED  FOR  END 
EFFECT  IN  LINE 

O — O THEORETICAL  KING- 
MIDDLETON  SECOND 
ORDER  IMPEDANCE  FOR 
CONSTANT  RADIUS 


* 


FIG.  3 


MEASURED  AND  THEORETICAL  IMPEDANCE  OF 
A HALF-DIPOLE 


FIG  4 COMPARISON  OF  EQUIPMENT  PROSES  WITH  DIPOLE  ANTENNA  THEORIES 


APPA  ?ENT  CAPACITANCE  AND  REACTANCE  OF  THE  GAP 
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Chapter  IV 

EXPERIMENTAL  DATA  AND  COMPUTATIONS 


1 Extent  of  Parameter  Variation 

The  parameters  selected  to  be  varied  m the  measurements  are  the 
lengths  of  the  driven  and  parasitic  elements  and  the  spar:  .r;.g  between 
them.  Since  one  of  the  purposes  for  making  the  measurements  is  to 
check  the  theory  over  its  complete  range  of  applicability  these  parameters 
must  be  varied  over  a range  which  includes  poor  as  well  as  good  agreement 
between  theory  and  experiment.  The  trial  functions  used  are  such  that 
the  current  distribution  is  well  represented  on  the  portion  of  the  antenna 
for  which  |z|  < g where  g is  less  than  \/Z.  Likewise  the  trial  current 
on  the  section  with  |z|  >g  is  a good  approximation  for  h g less  than  about 
3X./4.  Hence  half  lengths  ranging  up  to  a half  wavelength  for  the  driven 
antenna  and  overall  lengths  up  to  3X./4  for  the  parasite  should  be  investigated 
experimentally.  The  impedance  and  current  distribution  for  similar  lengths 
should  be  calculated  from  the  theory.  Longer  lengths  should  be  investigated 
experimentally,  but  there  is  little  point  to  computing  theory  for  them.  The 
theory  should  be  m good  agreement  with  measurements  for  short  lengths, 
g<  \ /4  and  h-g  less  than  \/2,  and  should  become  progressively  less 
applicable  at  greater  lengths.  The  theoretical  impedance  should  be  m 
better  agreement  with  the  measured  impedance  over  a greater  range  of 
lengths  than  that  over  which  the  theoretical  current  distribution  is  m 
agreement  with  the  experimental  distribution.  This  is  a consequence  of 
the  use  of  a variational  principle  to  improve  the  ’mpedance  approxima- 
tion. 

The  gap  should  be  varied  experimentally  from  zero  fi.e.  actual 
contact  of  the  elements)  to  as  large  a spacing  as  *.s  pcss  ble  in  the 
existing  experimental  equipment.  The  theory  however  cons-ders  the 
gap  to  be  well  represented  by  a variable  Jumped  capacitance  and  assumes 
the  coupling  between  the  cylindrical  surfaces  cf  the  array  elements  to 
be  constant.  Hence  the  limit  of  applicability  of  the  theory  is  at 
the  point  where  the  coupling  between  the  cylindrical  surfaces  has 
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change  ci  appreciably.  Note  also  that  the  coupling  between  the  end 
surfaces  changes  very  rapidly  with  spacing,  and  as  soon  as  this  change 
is  essentially  complete,  then  there  will  be  no  further  appreciable 
effects  due  to  actually  reducing  it  to  zero.  With  this  m mind,  the 
theory  is  limited  to  the  maximum  gap  size  required  to  reduce  the  end 
cap  coupling  effect  to  essentially  zero.  As  the  gap  is  increased 
beyond  either  of  the  above  limits,  whatever  effects  result  will  not  be 
predicted  by  theory. 

The  importance  of  these  limits  and  the  rapidity  with  which  they 
occur  are  difficult  to  predict.  A 1/10  wavelength  gap  should  make  a 
very  large  change  in  the  coupling  between  the  cylindrical  surfaces, 
so  the  theory  should  certainly  not  be  used  beyond  this  point.  Probably  i/20 
to  1/  50  wavelength  would  be  a reasonable  limit.  As  far  as  the  gap 
coupling  is  concerned,  consideration  of  Fig.  5 shows  that  the  gap 
capacitance  varies  only  very  slowly  for  gaps  greater  than  1 cm,  or  about 
1/50  of  a wavelength.  Beyond  this  point  the  gap  reactance  does  not 
change  appreciably  and  the  theory  is  no  longer  applicable.  Hence  the 
gaps  to  be  investigated  by  computation  from  the  theory,  and  experimentally, 
should  range  at  least  over  separations  of  zero  to  1 cm.  Larger  gaps  may 
be  investigated  experimentally  as  far  beyond  as  is  convenient. 

2.  Measured  Impedance  Data. 

The  impedances  measured  are  shown  in  Figs  6 through  IS  along 
with  the  corresponding  arrays  on  which  the  measurements  wets  made. 
Smith  charts  were  chosen,  in  Figs.  6 through  11  for  example,  rather 
than  rectangular  plots  since  the  curve  of  the  collinear  impedance  between 
the  two  end  points  would  be  more  obvious.  When  the  circle  was  too  small 
to  be  convenient  the  impedances  were  denormalized  and  plotted  in  rectan- 
gular resistance  and  reactance  coordinates, 

3 Computational  Procedure. 

The  choosing  of  cases  for  computation  must  be  done  with  care, 
for  the  calculations  are  complicated  and  time  consuming.  On  this  basis, 
it  is  convenient  to  consider  cases  for  which  as  many  terms  as  possible 
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are  zero  in  the  expressions  for  the  "y's  used  in  the  impedance  expres- 
sion. The  trigonometric  functions  contained  therein  are  zero  for  combina- 
tions of  h and  g equal  to  multiples  of  a quarter  wavelength  The.  other 
important  aspect  is  to  choose  a series  of  lengths  that  veil  check  the 
theory  over  at  least  its  expected  range  of  use.full.ne  s s . As  described 
in  sectionJV  1 such  lengths  are  for  g up  to  \ /2  and  h~g  up  to  3V  /4. 

Such  a series  of  combinations  are  listed  in  the  table  on  page  M In  addition 
are  listed  the  end  point  dipoles  into  which  the  array  degenerates  for  gaps 
equal  to  zero  and  infinity;  the  pertinent  figures  are  also  given.  Figure  39 
is  a plot  of  the  measured  half  dipole  impedances  on  the  line  with  a few  of 
the  array  impedance  spirals. 


The  current  distribution  for  a particular  configuration  is  computed 
using  Eq.  (2-3)  after  first  evaluating  (2  26)  and  {2  23)  for  Ij(z);  Ig(z)  is 
obtained  using  the  equations  on  p 24  The  phase  and  amplitude  of  the 
current  distribution  have  been  plotted  as  listed  on  page  14.  The  experimental 
data  have  been  plotted  by  adjusting  the  experimental  amplitude  at  the  peaks 
to  be  the.  same  as  for  the  theoretical  curves  with  the  exception  of  Fig.  41 
which  has  been  plotted  so  the  amplitude  at  z = n «-  or  respond  8 to  the  measured 
driving -point  impedance  for  this  case.  The  position  of  the  current  relative 
to  the  distance  scale  has  been  occasionally  adjusted  to  account  for  errors  in 
knowing  the  position  of  the  auxiliary  probe  on  the  polyfoam  column  supporting 
the  array. 


Equation  (2-4)  for  the  theoretical  array  impedance  requires  the 
use  of  the  current  distribution  at  z = 0 and  at  z = g.  The  impedance 
Z is  the  apparent  reactance  X of  the  gap.  may  be  computed  directly, 

but  there  is  occasionally  some  additional  information  which  may  be  either 
directly  incorporated  into  the  computations  or  merely  used  as  a check 
on  the  infinite  gap  The  value  of  Z^  for  the  infinite  gap  point  is  occasionally 
known  from  another  computation,  for  example,  for  configuration  4 of  the 
table  Z (gap  = oo)  is  that  for  a dipole  of  half  height  equal  to  X.  / 2.  This 
impedance  is  known  since  it  is  a necessary  part  of  the  computation  of 
Zc  for  configuration  2.  It  may  be  used  as  follows  Equation  (2-4)  is 
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Evaluating  this  for  the  experimentally  infinite  gap  point  or  for  Z = co 
in  (2-4)  there  results 


Zc(gap  = oo)  = Zs(h)/|[  1 - DiO$S(g)j  (4  1) 

However  this  Z^(gap  = oo)  is  that  of  a singly  driven  dipole  of  half  height 
equal  to  the  gap  position  g;  that  is. 


Z (gap  = oo)  = Z (half  height  = g)  (4-2) 

c s 

The  symbol  Z^half-height  = g)  is  used  to  avoid  any  conlusion  resulting 
from  using  the  letter  h in  a description  of  this  dipole.  Substituting  (4-2) 
in  (4-  1 ) yields 


D(Q)S(g)  = 1 


Z (h) 


(naif  height  = g) 


(4-3) 


Equation  (4-3)  was  used  in  the  computation  for  Fig.  17  and  considerably 
improved  the  theoretical  agreement  with  experiment  It  was  also  used  in 
Fig.  16  and  Fig.  18  but  had  very  little  effect  since  the  agreement  was 
already  quite  good.  It  was  not  used  in  Fig.  15.  Note  that  Z^  (half-height  = g) 
is  obtained  theoretically  by  evaluation  of  2-29  for  the  height  of  interest. 


4.  Conclusions 

The  agreement  between  theoretical  computations  and  experimental 
measurements  is  as  considered  in  Section  IV- 1 When  g<  (h-g)  and  both 
are  short,  the  measurements  are  in  best  agreement.  The  agreement  is 
generally  good  for  configuration  1 (see  table  on  page  34)  for  a . 2X. 
driving  element  and  a . 3\  parasite.  It  is  less  so  for  configuration 
2 in  which  both  elements  are  X.  / 4 long,  and  it  is  least  good  for  the  4th 
case.  The  impedances  generally  agree  quite  well  for  all  except  configura- 
tion 4 for  which  the  conditions  on  the  trial  current  distributions  are 
pot  met.  This  is  to  be  expected  since  the  variational  principle  technique 
allows  the  use  of  a rather  poor  current  approximation  in  obtaining 
comparatively  good  impedance  values.  Note  that  the  agreement  (see  Fig  18) 
of  the  theory  of  King  is  good  for  large  spacing s if  the  theor  0 tic  oil  curve  were 
displaced  such  that  the  infinite  gap  point  was  superimposed  onto  a better 
end  point  impedance  than  that  obtained  in  his  zero-order  theory.  The 
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discrepancy  for  small  spacings  is  to  be  expected  since,  his  theory  is  based 
entirely  on  varying  the  coupling  between  the  cylindrical  surfaces  of  the 
antenna  elements  and  will  be  in  error  when  the  end-cap  coupling  is  appreciable 

On  the  Smith  Chart  of  Fig.  39  is  plotted  a curve  for  the  measured  half 
dipole  impedance  and  a few  of  the  impedance  arcs  for  eollinear  arrays.  This 
makes  the  relation  of  the  end-point  impedance  of  the  array  to  those  of  the  half 
dipole  quite  obvious.  It  also  makes  it  possible  to  guess  roughly  where  an  im- 
pedance would  fall  for  some  other  configurations  than  those  considered  here. 

In  conclusion,  a theory  for  the  close-spaced  eollinear  array  has 
been  formulated  on  the  basis  that  the  array  is  the  superposition  of  a 
doubly  and  a singly  driven  dipole.  The  theory  has  been  evaluated  for  a 
series  of  specific  configurations  that  arc  typical  of  the  applicable 
range  of  the  theory.  These  same  configurations  have  been  investigated 
experimentally  as  a check  on  the  basic  assumptions  in  the  theory  and 
also  as  a check  that  the  region  of  usefulness  of  the  theory  has  been  properly 
estimated. 

5.  A Two  -Wire  Line  as  a Coupling  Reactance  between  Elements 

This  type  of  coupling  has  been  very  briefly  considered  experimentally. 
The  devices  pictured  in  Fig.  40  have  been  used  to  couple  the  ends  of  the 
driven  element  and  the  parasitic  elem1  nts.  These  permit  the  coupling 
reactance  between  the  elements  to  be  varied  without  changing  their  spacing. 

The  experimental  setup  is  similar  to  that  used  elsewhere  except  small  two- 
wire  line  is  used  at  the  gap.  It  has  a characteristic  impedance  of  135  ohms 
and  a spacing  of  0.  2 inches.  Two  arrangements  were  used.  One  of  these 
is  a line  whose  overall  length  is  constant  and  upon  which  a shorting  bar  is 
moved.  The  other  is  one  the  overall  length  of  which  is  varied.  Both  are 
needed  for  a complete  study  of  the  problem.  T a first  wiii  have  a constant 
length  for  the  unbalanced  mode  on  the  open- wire  line;  the  other  will  present 
a variable  length  to  the  unbalanced  mode.  Plots  are  given  in  the  last  figure 
showing  the  measured  driving -point  impedance  using  both  of  the  line  types. 

The  points  labeled  constant  overall-length  line  are  those  for  the  line  using 
the  shorting  bar„ 


A= 50  cm 


FIG  40  MEASURED  IMPEDANCE  USING  TWO-WIRE  LINE 
COUPLING  BETWEEN  ELEMENTS 
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Appendix  A 

EVALUATION  OF  THE  y-INTEGRALS 


The  substitution  of  (2-26)  and  (2-27)  into  (2-25)  yields* 


■>£ 


Z - 
d *±ir 


-jw(l+aC)  / dz[  1 + C(a  + cos  (3z)]  / [ K(z,z')  + K(z,-z')]  dz' 


dz[  1 + C(a  + cos  (3z)j 
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[K(z.z')  + K(z,-z')]L(z')I  (z’)dz' 


Ah 


-ju(6  + €D)/  dz  [1  + C(a  + cos  pz)J 


[K(z,z*}  -!•  K(z,-z')]  ds 


-juj(  1 + aC)  / dz 
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-)] 


l-cosp(h-z)^  + D(  sin  p(n-z)  + £ 


1 -C03  (S( 
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[ K(z,z')  + K(z, 


z')]  dz' 


•>h 

+ j dz 
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6 |l -cos  3{h-z)j  + D(sinp(h-z)  + £ 1 1 -cos  j3(h- z) 
g+ 

f [K(z.z')  + K(z,-z')]  L(z-)Id(z-)dz' 
Jo  - 


::Most  of  t.he  following  expressions  are  taken  from  reference  6 
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dz[  6 ^'1  -cos  p(h-x)j  + D( sin p(h- z)  + * f 1 -cos  p(h- z)j  )] 


[ K(z,z’)  + K(z,-z')]dz' 


Jg 


6 

f°f 


[ l + C(a  + cospz)]  [K(z,h)  + K(z,-h)]  dz 


h 

- jjp-  D [ sjl-cos  p(h-  z)J  + D(sin|J(h-z)  + t ^ 1-cos  p(h- s)|  )] 

•) 

[K(z,h)  + K(z,-h)]  dzj  . 

Note  that  the  integral  in  the  variable  z*  may  be  simplified  by 
f+ 

[K(z.z')  + K(z,-z!)]  L(z')Id(z')dz*  = 
g- 


( A- 1 ) 


O - 

- -jy  [K(z,g)  + K(Zi-g)]  f (Id(z'))dz'  = 


p 


8+  2 

r 

if-  dz,i 


^r[  K(z,g)  + K(z,-g)]  [ e + D + C sin  fig] 
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Using  this  in  (A-l)  yields 


where 


Zd  =iw'  +YCC  + YDD  + YCDCD  +V  CCC  +7DDD  ^ 
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r (a  + cos  pz)[K{z,g)  + X(z,~g)]  dz 


g r 

+ ep  dz(o  +cospz)  j [K(z,z'}  + K(z,-z')]  dz1 


+ a (3  / 

v-4 y 


h g 

dzsinp(h-z)  C [K(z,z')  + K(z,-z')]  dz' 


i 


r 8 

+ at  p j dz  |l -cos  p(b  - z)j  C [K(z,z')  , K(z,-z')]  dz' 
Jo  Jr> 


+ sinpg  / sinp(h-z)[  K(z,g)  + K(z,-g)]  dz 


+ € sin  pg 


{ 1 -cos  p{h  z)}  [ K(  z,g)  + K(  z, -g)]  dz 


g 

t (a  + cospz)[  K{ z,h)  + K(z,-h)]  dz 


dz(a  + cospz)  /[  K(z,z')  + K(z,-z')]  dz' 


+ sin  Pg  / (a  + cos  pz)[  K(z,~  g)  + K(z:g)]  dz 

Jo 

g h 

[ K(z,g)  + K(z,-g)]  dz  + J'  [K(z,h)  + K(z,-h)]  dz 

h g 

+ tp  f dz  | 2-cos  p(h- z)|  / [K(z,z!)  + K(z,-z‘)J  dz' 


v'o 


\ 


D 


TR178 


-41- 


g 

. r r 


+ P j dzsinp{h -z)  j [ K(z,z’)  + K(z,-z')]  dz' 

Vn 


g 


h 

+ e f sin  |3{h- z)[  K(  z,g)  + K(z,-g)]  dz 


T ( 6 + € 


{1-cos  p(h-z)]  [K(zlS>  + K(z,-g)]  dz 


n h 

+ 5?  f dzsinjj(h-z)  j [K(z,z!)  + K(z,-z’)]  dz' 

Jg  Jg 


h h 

+ 26c  £ / dz  { 1 -cos  p(h- z)\  C [K(  z,z')  + K(z,-z')l 

Jg  ' 


J I 


DD 


h 

+ 6 / {l-cos  p(h-z)|  [ K(z,h)  + K(z,-h)]  dz 
h 

= f sinp(h-z)[K(z;g)  + K(z,-g)]  dz 


+ c 


(l-cos  P(h  zjjlKtz^g)  r K(z,-g)]  dz 


h h 

+ «p  C dzsinp(h-z)  f [K(z,z')  *■  K(  z , - ) j dz: 


X 


h h 

■5-  c f dz.{l -cos  p(h-z)j  j [ K(  z,z')  + K(z,-z')]  <tz* 
Jg  Je 


+ J ain*ph(h-z)^K(2.,h)  + K{z, -hg  d* 


jg 


dz' 


TR178 


-42- 


h 

+ £ / |>-cosp(h-z)]  [K(z,h)  + K(z,-h>]  dz 


The  above  y expressions  contain  integrals  of  six  different  forms. 
Three  may  be  evaluated  as  follows: 
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an  integration  by  parts  yields 
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Another  form  is 
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Ji  Ji 


c 


c 


-f 


cos  pz  K(  z,g)  dz 


cos  pz  K(z,£)  dz 


G(b,c,f,g)  = cos  pb[  C(pa,p(g-b)  - C(pa,p(f ~b))] 

- cos  pc[Ctpa,p{g  -c)  - C{pa^p(f-c))j 

- cos  pg[  Cc(pa,p(c-g)  - <^c{pa,p(b-g))] 


+ •inpg[C8(pa,p{c-g))  - Cypa.pic-g))] 
+ cos  Pi^  ^^{pa,p(c -f))  = ^c(?a,p{c-f))] 

- »in pf[  Cg(pa,p(c-f))  - Cs(pa,p{b-.f))] 

-j  f cos  pb[  S(pa,P(g  -b))  - S<pa,p{f-b))] 

- cospc[S{pa,p(g-c))  - S(patp(f-c))] 

- cos  pgf  Sr{pa,p{c-g))  - Sc(pa,p(b-g))] 


+ »inpg[Sa(pa,p(c-g))  - Sg(pa,p(b-g»] 


I V.UO  j 


>spf[Sc(Pa.6{c-f»  - Sc(pa,p(b-f))3 
- sin  pf[  S (pa,p(c -f))  - 3 <«a„3(b-f»]  ] . 

3 3 J 

By  a similar  process 


c 

r> 


H(b,c,f,g)  = 0 j dzcospz  / K(z,z')  dz 

Jb  Ji 
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H(b,c,f,g>  = sinpc[  C(pa,p(g-c))  - (?(pa,p(f-c))] 

- flinpb[C(pa,p(g-b»  - C(pa,p(f-b))] 

+ cos  Pg{  Gs(pa.p(c-g))  - Cg(pa,p(b-g;)j 

+ sinpg[  Cc(pa,p^c-fe))  - C (pa.p(b-g))] 

- cospf[Cg(pa,p(c-f))  - Cg([pa,p\b-f))] 

• 3inpf[Cc(pa,p(c-f))  - Cc(pa,p(b-f»] 

-j  [ sinpc[  S(pa,p(g-c))  - S(pa,p(f-c))] 

- smpb[S(pa,p(g-b))  - S(pa,p(f-b))] 

+ cos  pg[  Ss(pa,p(c-g))  - Ss(pa,B(b-g»] 

+ sinpg[  Sc(pa,p(c-g))  - Sc(pa,p(b-g))j 

- cospf[S  (pa.p(c-f))  - S (pa.p(b-f))] 

8 8 ij 

- sinpf[  Sc{pa,p(c-f»  - Sc(pa,p(b-f»]  } 

Three  single  integrals  are  immediately  expressable  in  terms  of  the 
tabulated  functions.  These  integrals  are 

c 

I(b,c,g)  = / cos  pz  K(-s,g)  dz 

v4> 

= cospg[  Cc(pa,p(c-g))  - Cc(pa,p(b-g))] 

- sinpg[  Cg(pa,p(c-g))  - Cg(pa,p(b-g))] 

-j  ( cospg[  Sc(pa,p(c-g))  - Sc(pa,p(b-g))] 

- sinpg[Sg(pa,p(c-g))  - Sg(pa,p(b-g))]  J 
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sinpg  K(z,g)  ds 


= cos  pg[  Cs(pa,p{c-g^  - Cs(8a,p(b-g))] 
+ sin?g[  Cc(pa,§(c-g»  - Cc(pa,8fb-g|D3 
-j  ( cospg[SsCpa,pic-g))  - Saipasp(b-g])D| 

+ sinpg[  Sc(pa,p(c-g))  - ScfRa,p(b-gHl  j 
c 

M(b,c,g)  = C K(z,g)  dz 

J b 


= <T(pa,p(c-g))  ^(pa.p(b-g)) 

- j[s(Pa,p(c-g))  - S(pa,p(b-g))| 

The  following  symmetries  are  applicable  throughout 


P 


g 


K(z,-z')dz' 


-F(b,c,-f,-g^ 


P j dz  sinpz  / K(z,-a')dz'  = -G(b„c ,-f,-g) 

Jb  Jl 


g 

dz  cos  p z f 

Jf. 


K(z,-z'}dz'  = - H(b,c,-f,-g) 


The  functions  S,  S , S , C , £T  are  defined  and  tabulated  very  completely 

S C 8 C 

in  the  reference.  All  are  odd  about  the  origin  except  S and  C - 

s s 

Note  that  in  evaluating  the  various  integrals  that  there  is  a 
definite  relationship  between  the  tabulated  functions  appearing  in  real 
parts  of  F,G,  etc.  and  those  that  appe^^  in  the  imaginary  parts  of  the 
integrals  One  of  these  for  example,  is 


H = f [C  ( wj,  (J  (x).  C(v>,  sinz]  - jf[S  (w),  S lx),  S{y),  -cor  z] 

SC  *sc 


The  real  parts  of  the  integrals  contain  only  Cs«  Cc>  C,  and  sine  functions 

while  the  imaginary  parts  contain  an  identical  arrangement  of  S , S , S, 

s c 

and  minus  cosine  functions. 
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